Dynamical spin susceptibility and the resonance peak in the pseudo-gap region of the 

underdoped cuprate superconductors 
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We present a study of the dynamical spin susceptibility in the pseudogap region of the high-T c 
cuprate superconductors. We analyze and compare the formation of the so-called resonance peak, 
in three different ordered states: the d a .2_ H 2-wave superconducting (DSC) phase, the d-density wave 
(DDW) state, and a phase with coexisting DDW and DSC order. An analysis of the resonance's 
frequency and momentum dependence in all three states reveals significant differences between them. 
In particular, in the DDW state, we find that a nearly dispersionless resonance excitation exists only 
in a narrow region around Q = (n, n). At the same time, in the coexisting DDW and DSC state, the 
dispersion of the resonance peak near Q is significantly changed from that in the pure DSC state. 
Away from (n, ir), however, we find that the form and dispersion of the resonance excitation in the 
coexisting DDW and DSC state and pure DSC state are quite similar. Our results demonstrate 
that a detailed experimental measurement of the resonance's dispersion allows one to distinguish 
between the underlying phases - a DDW state, a DSC state, or a coexisting DDW and DSC state - 
in which the resonance peak emerges. 

PACS numbers: 71.10.Ca,74.20.Fg,74.25.Ha,74.72.-h 



I. INTRODUCTION 

One of the most controversial topics in the field of 
high-temperature superconductivity is the origin of the 
so-called 'pseudogap' phenomenon observed by various 
experimental techniques in the underdoped cuprates (for 
a review see Refi and references therein). A large num- 
ber of theoretical scenarios have been proposed to ex- 
plain the origin of the pseudogap 2 * 3 ". Among these is the 
d-density wave (DDW) scenario 3 - which was suggested to 
explain some of the salient features of the underdoped 
cuprates such as the d a .2_ y 2-wave symmetry of the pseu- 
dogap above T c , the anomalous behavior of the superfiuid 
density 4 and of the Hall number—, as well as the presence 
of weak (orbital) antiferromagnetism 6 . The DDW-phase 
is characterized by circulating bond currents which alter- 
nate in space, break time-reversal symmetry and result in 
an orbital (antiferromagnetically ordered) magnetic mo- 
ment. 

In this article, we investigate the momentum and fre- 
quency dependence of the dynamical spin susceptibil- 
ity x(q, w), in the underdoped region of the cuprate 
superconductors. In particular, we compare the forma- 
tion of a resonant spin excitation (the "resonance peak" ) 
in three different ordered states: the DDW phase, the 
d x 2_ y 2-w&ve superconducting (DSC) phase, and a phase 
with coexisting DDW and DSC order. The observa- 
tion of the resonance peak in inelastic neutron scatter- 
ing (INS) experimenta 7 i 8 i 9 i 10 i n i 12 i 13 i 14 i 15 is one of the key 
experimental facts in the phenomenology of the high- 
Tc cuprates. In the optimally and overdoped cuprates, 
the resonance peak appears below T c in the dynamical 



spin susceptibility at the antiferromagnetic wave vector 
Q = (77,71"). In the optimally doped cuprates, the res- 
onance's frequency is u> res « 41 me\2i 1 i, a frequency 
which decreases with increasing underdopingi 2 * 1 ^. A 
number of theoretical scenarios have been suggested for 
the appearance of the resonance peak in a superconduct- 
ing state with d x 2_ y 2-wa.ve symmetryi£iiLi£. In one of 
them, the so-called 'spin exciton' scenarioiSii 7 ., the reso- 
nance peak is attributed to the formation of a particle- 
hole bound state below the spin gap (a spin exciton), 
which is made possible by the specific momentum depen- 
dence of the d x 2_ y 2-wave gap. Within this scenario, the 
structure of spin excitations in the superconducting state 
as a function of momentum and frequency can be directly 
related to the topology of the Fermi surface and the phase 
of the superconducting order parameter. This scenario 
agrees well with the experimental data in the supercon- 
ducting state of the optimally and overdoped cuprates. 
In the underdoped cuprates the resonance-like peak has 
also been observed in the pseudogap region above T c as 
well as in the superconducting statsiS* 1 ^* 1 ^ 5 -. In this ar- 
ticle, we address the question whether in the underdoped 
cuprates, the resonance peak above T c emerges from the 
presence of a DDW state, as first suggested by Tewari et 
alX , and below T c from the coexistence of a DSC and 
DDW phase. To answer this question, we develop a spin 
exciton scenario for the pure DDW-phase as well as the 
coexisting DSC and DDW phases. By studying the de- 
tailed momentum and frequency dependence of the res- 
onance peak in both phases, and by comparing it with 
that in the pure DSC state, we identify several charac- 
teristic features of the resonance peak that allow one to 
distinguish between the underlying phases, in which the 
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resonance peak emerges. 

The remainder of the paper is organized as follows: 
in Sees. ITT1 and IIIII we discuss the form of the resonance 
peak in the pure DDW phase and the coexisting DDW 
and DSC phase, respectively, and compare it with that in 
the pure DSC state. In Sec. lIVI we summarize our results 
and conclusions. 



II. PURE DDW STATE 

Starting point for our calculations in the pure DDW- 
state is the effective mean field Hamiltonian 



H 



J 'OW — ^ £kC kiCT Ck,er + ^ 
k,<r k,er 



iM/ k C k l( r C k+Q^ I 1 ) 



where Q = (it, it) is the ordering wavevector of the DDW 
state, Wk — ^-(coskx — cosk y ) is the DDW order pa- 
rameter, 



ck = — 2i (cos k x 



cos k y ) — At' cos k x cos k y 



■H (2) 



is the normal state tight-binding energy dispersion with 
t, t' being the hopping elements between nearest and 
next-nearest neighbors, respectively, and \i is the chem- 
ical potential. In the following we use t = 250meV, 
t'/t = -0.4 and /z = -1.083i. The Fermi surface (FS) 
obtained from Eq. describes well the FS measured by 
photoemission experiments on Bi2Sr2CaCu20s+^. In 
order to directly compare the dynamical spin suscepti- 
bility in the DDW state with that in the DSC state, we 
take the DDW order parameter, Wo = 42meV, to be 
equal to that in the DSC stated. We note here, that the 
above Hamiltonian can be obtained from a microscopic 
Hamiltonian with short-range repulsion or superexchangc 
interactions22i2i. After diagonalizing the Hamiltonian, 
Eq.JTJ, on e finds that the excitation spectrum possesses 
two bands with energy dispersion 



E, 
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(3) 



where = (ek±£k+Q)/2. In Fig^a) we present the re- 
sulting Fermi surface in the DDW phase. Due to the dou- 
bling of the unit cell in the DDW state, the Fermi surface 
consists of hole pockets centered around (±7r/2, ±7r/2) 
and electron pockets around (±7r,0) and (0, ±7r). This 
type of Fermi surface has not yet been observed exper- 
imentally in the underdoped cuprates, possibly, as has 
recently been argued, due to additional interactions be- 
tween quasiparticles^. For the above band parameters, 
the chemical potential lies within both branches of the ex- 
citation spectrum thus preventing the formation of a gap 
at the Fermi level. This is clearly visible from Fig^b) 
where we plot the density of states (DOS) for various 
values of the DDW gap. In particular, one finds that 
for \t'\ > Wo/A a suppression of (i.e., dip in) the DOS is 
formed away from the Fermi level which increases with 
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FIG. 1: (color online) (a) Fermi surface in the DDW state for 
t = 250meV, t'/t = -0.4, fj, = -1.083* and Wo = 42meV. 
The black and gray arrows represent the magnetic scattering 
vectors for Q = (n, n) and Qi = 0.98(-7r, it), respectively, (b) 
Calculated DOS for several values of the DDW gap. 



increasing W - In contrast, for \t'\ < Wo/4, this sup- 
pression, which we identify with the pseudogap, opens at 
the Fermi level as was noted previously^. Note that for 
t' = 0, the DOS vanishes at the Fermi level, and the DOS 
resembles that of a d x 2_ y 2-weive superconductor. 

In order to compute the dynamical spin susceptibility 
in the DDW state, we first introduce the spinor 



k.tT 



L k,(T' °k4 



Q,(T 



(4) 



where a is the spin index, and the electronic Greens func- 
tion in the DDW state is defined as G^k, r — r') = 
— (T\E , k,cr('r) x I'j <; ff (r')). The bare (non-interacting) part 
of the dynamical spin susceptibility per spin degree of 
freedom is then given by 



Xo(q,«^m) 



E' Tr [G (k ''' 



~ 8 



x G(k + q, iuj n — ifl,, 



(5) 



where G(k, iu> n ) — G a (/k,iuj n )ao is the Green's function 
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matrix in momentum and Matsubara spaced, and the 
primed sum runs over the reduced fcrmionic Brillouin 
zone of the DDW state. Note that with the above defi- 
nition, one has x§ z — 2xo- After performing the summa- 
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tion over the internal Matsubara frequencies and analytic 
continuation to the real frequency axis, one obtains for 
the retarded spin susceptibility in the DDW phase 



Xo(q,w) 
/ 

+ 1 

V 



-E' 



1 + 



g k £ k+q + WkW k +q 
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u + £+ +q + E+ ^ u> + z0+ - £ k+q + 

/(*£+,) -/TO \ 



W + »0+-2£ +q + i£ 
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where /(e) is the Fermi function. 

We first analyze the behavior of the imaginary part of 
Xo at Q = (-7T, it), and present in Fig. Gllnixo(Q, to) as a 
function of frequency in the normal state, the DSC state, 
and the DDW stated (for the form of xo in the DSC 
state, see Reftf). The behavior of Imxo(Q,^) in the 
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2: (color online) Imaginary part of xo(Q,<^) at Q = 
(n,n) in the normal, DSC, and DDW state as a function of 
frequency. The momentum dependence of the superconduct- 
ing gap is taken to be A(k) = Ao(cosfc I — cosfc y )/2, with 
Ao = Wo = 42meV. Inset: Imxo(Q, w) in the DDW phase as 
a function of frequency around Q® r DW . 

normal and DSC state have been extensively discussed 
in the literature (see, for example, Ref ^LiLiS) . In the 
normal state Imxo increases linearly at low frequencies 
with a slope determined by the Landau damping rate, 
while, at higher energies its behavior is determined by 
the presence of the van Hove singularity. In contrast, 
in the superconducting state the susceptibility is gapped 
up to an energy ^ r sc — mink(|Ak| + |Ak+cj|) where 
Ak is the superconducting gap and both k and k + Q 
lie on the Fermi surface. Due to the symmetry of the 
superconducting gap, one finds Ak = — Ak+Q, resulting 
in a discontinuous jump of Imxo at Q® r s ci6 i 17 . 

In order to discuss the behavior of Imxo in the DDW- 



state, we first note that the expression for xo in Eq.© 
contains two terms that describe intraband scattering 
within the .E^-bands, and two terms that represent inter- 



band scattering between the two bands. Since E. 



k+Q 



the intraband scattering terms do not contribute to 
Imxo at Q. Moreover, since E^ < E£ the first interband 
scattering term yields a non-zero contribution to Imxo 
only for negative frequencies. Thus, only the second in- 
terband scattering term in Eq. © contributes to Imxo at 
Q. Since the Fermi surfaces of the two energy bands, E^ , 
cannot be connected by the wave vector Q, as is evident 
from Fig^a), Imxo is gapped at low frequencies up to 
an energy Q^, w 64.8 meV. The latter is determined 

by the minimum value of E£ — E^ = 2 \J (e^) 2 + W£ in 
the DDW Brillouin zone, a condition that is set by the 
(^-function arising from the last term in Eq.© for Imxo- 
Note that e k = along the boundary of the DDW Bril- 
louin zone. Due to the requirement sgn(E^) ^ sgn(E^), 
we find that il® r DW = 2| Wk 1 where ko is the momentum 
at which the hole pocket around (tt/2, it/2) is intersected 
by the DDW Brillouin zone boundary (see Fig.^l. More- 
over, since £k+Q = — e k an< ^ ^k+Q = — Wk, the second 
coherence factor in Eq.© is identical to 2 for all mo- 
menta. Note that there exist two important differences 
in Imxo between the DDW and DSC state. First, in 
the DSC state, Imxo 7^ requires that the frequency ex- 
ceeds Qcr SC = mink (l^k| + |Ak+Q|), a condition which 
is set by the <5-function in Imxo (see Eq.©) and simply 
reflects energy conservation. In contrast, in the DDW- 
state, Imxo ^ requires (a) that lj — ^k+Q + = ^ f° r 
certain momenta k, and (b) that for the same momenta 
/(^k+q) - /C^k ) 7^ 0. We find that there exist momenta 
for which (a) is satisfied at frequencies u < n° DW , but 
that for the same momenta f(E£ + ) — f(E^) = (at 

T = 0). In other words, the critical frequency Q® r DW 
for the onset of a non-zero Imxo is determined by the 
difference in the population of the states that are in- 
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volved in the scattering process, and not by energy con- 
servation as in the superconducting state. This quali- 
tative difference between the DDW and the DSC state 
bears important consequences: for uj > il® r DW , one has 
Imxo ~ x/^ — fl®. DW , in contrast to the discontinuous 
jump of Imxo at fl® r sc in the DSC state (this result 
for the DDW state differs from that in Refi^ due to the 
different Fermi surface topology considered here). This 
behavior becomes immediately apparent when one plots 
Imxo in the DDW-state around ui — fl^ r DW , as shown in 
the inset of Fig. [21 Note that the number of momenta 
which are involved in scattering processes and thus con- 
tribute to Imxo rapidly increases for ui > £l® r DW due 
to a steeply increasing density of states of the function 
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FIG. 3: (color online) Imxo in the DDW state as a function of 
frequency for several momenta q = r](n,Tt) along the diagonal 
of the magnetic Brillouin zone (BZ). Inset: Imxo in the DSC 
state as a function of frequency for several momenta q — 
77(71", tt). 



For momenta q ^ Q, the behavior of Imxo in the 
DDW state is more complex, since in addition to inter- 
band scattering, intraband scattering is now possible. In 
Fig. |3 we plot the frequency dependence of Im^o for sev- 
eral momenta in the DDW state (for comparison, Imxo 
in the DSC state is shown in the inset). We find that 
as one moves away from Q = {it, 7t), several square- root- 
like increases in Im%o appear, with the one lowest in en- 
ergy rapidly decreasing in frequency. In order to better 
understand the combined frequency and momentum de- 
pendence of Imxo m the DDW state, we present in Fig.0] 
the contributions to Im^o at Q; = 0.98Q from interband 
scattering [Fig.^a)] and intraband scattering within the 
.E^-band [Fig.^b)] and -band [Fig.QJc)] separately. 
Note that while the contribution from intraband scat- 
tering is approximately two orders of magnitude smaller 
than that from interband scattering, the former contin- 
uously increases from zero energy, such that Im%o does 
not any longer exhibit a gap. This result is valid for 
all momenta q ^ Q in the vicinity of Q. In contrast, 
the interband scattering term possesses three critical fre- 



0.20 
0.15 
0.10 
0.05 - 



0.0000 



(a) 

interband scattering 



(3) 



(2) 




intraband scattering 
Ek-band 



intraband scattering 
Ek-band 



10 20 30 40 50 60 70 80 90 100 

Frequency [meV] 



FIG. 4: Contributions to Imxo at Q 4 = 0.98Q from (a) inter- 
band scattering, (b) interband scattering within the iJ^-band, 
and (c) interband scattering within the S^-band. The arrows 
in (a) indicate the critical frequencies for the opening of in- 
terband scattering channels. The transitions in the fermionic 
BZ corresponding to the opening of the interband scattering 
channels are shown in Fig^a). 



quencies, Q, { ${i = 1,2,3), which arise from the opening 
of three non-degenerate scattering channels that are de- 
scribed by the scattering momenta shown in Fig. ^a). 
These three critical energies are indicated by arrows in 
Fig. 0{a). For all three scattering channels, the coher- 
ence factor is approximately 2. Note that the first and 
third scattering channel, which open at « 43 meV 
and fiir ~ 88 meV and are described by arrow (1) and 
(3) in Fig. ^a), respectively, connect momenta k and 
k' with k — k' = Qi — (tt, tt) and thus represent umk- 
lapp scattering. In contrast, channel (2) which opens at 

~ 76 meV [see arrow (2) in Fig. Qfa)] describes di- 
rect scattering with k — k' = Q; . The opening of each 
of these three scattering channels is accompanied by a 
square-root like increase of Imxo- Note that the lowest 
threshold frequency for interband transitions vanishes at 
Q,; = 0.94Q, since this wave vector connects momenta 
on the Fermi surfaces of the E£ and E^ bands. 

The emergence of a resonance peak in the DDW state 
can be understood by considering the dynamical spin 
susceptibility within the random phase approximation 
(RPA). Within this approximation^, the susceptibility 
(per spin degree of freedom) is given by 



x«PA(q,w) 



xo(q,w) 



(7) 



i - c^xo(q,^) 

where U is the fermionic four-point vertex. We first con- 
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sider q = Q and note that in the superconducting state, 
the discontinuous jump in Imx leads to logarithmic sin- 
gularity in Rexo- As a result, the resonance conditions, 
?7Rexo(Q,^ = w res ) = 1 and Imxo(Q,^ = w res ) = 0, 
can be fulfilled simultaneously below the particle-hole 
continuum for an arbitrarily small value of U > 0, leading 



to the emergence of a resonance peak as a spin exciton. 
In contrast, in the DDW state, Imx exhibits a square- 
root like frequency dependence above tt® DW , leading to 
an increase of Rexo a t the critical frequency, but not to 
a singularity in Rexo- Specifically, we find 



Rexo = 




where W = E c - Q® r DW , E c is the high energy cut-off 
for the square-root like frequency dependence of Im^o = 
ay/tu - n° DW , A± = tt° DW ±Lu and the ellipsis denote 
background contributions to Rexo that are independent 
of the opening of a new scattering channel. The above 
form of Rexo implies that U now has to exceed a critical 
value, U c , before a resonance peak (in the form of a spin 
exciton) can emerge in the DDW state. We note, how- 
ever, that the values of U typically taken to describe the 
emergence of a resonance peak in the DSC state of op- 
timally doped cuprate superconductors, exceed U c , such 
that a resonance peak also emerges in the DDW state. 
In other words, for U > U c , the resonance conditions 
[/Rexo(Q, w = w res ) = 1 and Imxo(Q, w = w res ) = are 
satisfied in the DDW state at a frequency uj res < fl%P w . 
As a result, a resonance peak emerges in the RPA spin 
susceptibility around Q = (7r,7r), as shown in Fig|5j (the 
value of U is chosen such that in the DSC state, uj res = 41 
meV). Away from Q, the mode becomes rapidly damped 
due to the opening of a scattering channel for intraband 
transitions, as discussed above. In addition, the low- 
est critical frequency for interband transitions rapidly 
decreases to zero. As a result, the resonance peak in 
the DDW state is confined to the immediate vicinity of 
Q = (tt, tt), and shows no significant dispersion. Note, 
that the upward and downward structures in ImxRPA 
visible in Fig. [S] do not represent real poles in the suscep- 
tibility but arise from the frequency structure of Imxo 
away from (tt,tt). This momentum dependence of the 
'spin exciton' in the DDW state stands in stark contrast 
to the dispersion of the resonance peak in the DSC stated 
(see Fig.EJ. 



III. COEXISTING DDW AND DSC PHASES 



We next consider a state with coexisting DDW and 
DSC order whose mean-field Hamiltonian is given by 



H 



k 



(9) 




FIG. 5: (color online) Imaginary part of the RPA spin sus- 
ceptibility in the DDW state as a function of frequency and 
momentum along the diagonal of the first BZ for U = 2.24eV. 



where ^ = (dp 4 +QT , c_ k-L , c_ k _ Q ^ and 4 



H k = 



e k iWu A k 

-iW k e k+Q -A k 

A k -e k iW k 

-A k -iWk -£ k +Q 



(10) 



The energy bands arising from diagonalizing the Hamil- 
tonian in Eq. i|10|) are given by 



A k 



(11) 



with e£ being the energy bands of the pure DDW state 
given above. The bare susceptibility, xo in the coexisting 
phase can again be calculated using Eq. © with the only 
difference that the Green's function G CT (k, icu n ) is now a 
(4 x 4) matrix. The full expression for xo in the coexis- 
tence phase is somewhat lengthy and therefore given in 
Appendix 

In Fig. we present Im^o as a function of frequency 
for several momenta q = 77(71", tt) along the diagonal of 
the magnetic BZ. At Q = (it, tt) (77 = 1.0), Imxo ex- 
hibits a single discontinuous jumps at the critical fre- 
quency, Sl c c ° ex = 97 meV. The magnetic scattering asso- 
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FIG. 6: (color online) Im xo as a function of frequency in 
the coexisting DDW and DSC state for several momentum 
q = 77(71", tt) along the diagonal of the magnetic BZ. 



ciated with the opening of this scattering channel con- 
nects the "hot spots" in the fermionic BZ, i.e., those 
momenta k and k + Q for which £k = £k+Q = 0. Cor- 
respondingly, the critical frequency is given by i\ c ° r &x = 
2y/ A 2 (ki ls ) + W 2 (ki ls ) where k; ls is the momentum of 
the hot spots. In contrast, away from Q = (tt, tt), we 
find that Imx exhibits 5 discontinuous jumps at critical 
frequencies, f2cr with i = 1,..,5 indicating the opening 
of new scattering channels (for r\ = 0.95 these five dis- 
continuous jumps are labeled in Fig. ■ Note that in the 
coexistence phase, the opening of a new scattering chan- 
nel is accompanied by a discontinuous jump, similar to 
the pure DSC state, but in contrast to the DDW state, 
as discussed above. The momentum dependence of these 

critical frequencies is shown in Fig. E[a). At Ocr [flcr ], 
a scattering channel for intraband scattering within the 
(^k ) band opens, and Im^o acquires a non-zero con- 
tribution from x (2) (x (5) ) given in Eq. (|X3|l [Eq.lJXBJl] of 
Appendix 1X1 As follows directly from Eas. l|A3|) and l|A6() . 
the coherence factors associated with these two scatter- 
ing processes vanish identically at Q = (tt, tt), and hence, 

no value for ttc r ' can be defined at this momentum. 
However, away from Q = (tt,tt), the coherence factors 
are not longer zero, and two discontinuous jumps ap- 
pear in Imxo that are associated with the opening of two 
new scattering channels at fl C r ■ Note that the magni- 
tude of the jumps at flir' 2 ^ increases as one moves away 
from Q = (tt, tt), which is a direct consequence of the 
increasing coherence factors. Similar to the supercon- 
ducting state the lowest critical frequency, ficr\ reaches 
zero at 77 = 0.8(tt,tt). In contrast, at f2cr' 4 ' 5 ^ new scat- 
tering channels for interband scattering between the 
and bands are opened. These three critical frequen- 
cies are degenerate at Q = (tt, tt), but this degeneracy is 
lifted for q 7- Q, as follows immediately from Fig. |7fa). 
For T] — 0.98, we present in Fig. |7fb) the scattering mo- 
menta that are associated with the opening of the above 
discussed five scattering channels. For completeness, we 
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FIG. 7: (color online) (a) Momentum dependence of the crit- 



ical frequencies Qcr (i = 1, .., 5). (b) FS of the pure DDW 
state and magnetic scattering vectors of the 5 scattering chan- 



nels that open at £1 



that the scattering vectors (1) and (2) describe intra- 
band scattering within the SI J and bands, while the 
scattering vectors (3), (4) and (5) represent interband 
scattering. The scattering vectors (3), (4), and (5) are 
identical to those present in the pure DDW state [for 
comparison, see Fig.l (a)]. 

In Fig. |H| we present the RPA susceptibility in the co- 
existing DDW and DSC phase. Similarly to the pure 
DSC state, the discontinuous jump in Imxo in the co- 
existence phase is accompanied by a logarithmic diver- 
gence in Rexo> which in turn gives rise to a resonance 
peak below the particle-hole continuum for an arbitrary 
small fermionic interaction. A comparison with the RPA 
susceptibility in the pure DSC state shown in Fig. Efb) 
reveals that the frequency position of the resonances peak 
varies more quickly with momentum (near Q = (tt, tt)) 
in the coexistence phase than in the pure DSC state. 
This difference becomes particularly evident when one 
plots the dispersion of the resonance peak both in the 
coexistence phase and the pure DSC state, as shown in 
Fig. We find that the difference in the dispersion of 
the resonance peaks is particularly pronounced around 
Q (0.95 < 77 < 1.05), with a more cusp-like dispersion 



also present the Fermi surfaces for the bands. Note, in the coexistence phase. This cusp follows the form of 
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FIG. 8: (color online) (a) RPA spin susceptibility as a func- 
tion of frequency in the coexisting DDW and DSC state for 
several momenta, (b) RPA spin susceptibility as a function 
of frequency in the pure DSC state for several momenta. 
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FIG. 9: (color online) Dispersion of the resonance peak in the 
coexisting DDW and DSC state, as well as in the pure DSC 
state. 



DSC state and the pure DSC state. However, away from 
Q, the particle-hole continuum, as well as the dispersion 
of the resonance peak are quite similar in both phases. 
IV. CONCLUSIONS 



In conclusion, we have analyzed the momentum and 
frequency dependence of the dynamical spin susceptibil- 
ity in the pure DDW state and the phase with coexisting 
DDW and DSC order. We find that due to the opening 
of a spin gap in Imxo in both phases, a resonance peak 
emerges below the particle-hole continuum. However, in 
the DDW state, Im%o exhibits a square-root like increase 
at the critical frequencies, in contrast to the coexisting 
DDW and DSC phases (or the pure DSC phase), where 
the onset of Im^o ^ is accompanied by a discontinuous 
jump. As a result, Rexo in the DDW state does not 
exhibit a divergence, but simply an enhancement at 
the critical frequency, and hence, a finite fermionic 
interaction strength is necessary for the emergence of a 
resonance peak in the DDW state. This result is qual- 
itatively different from the coexisting DDW and DSC 
phase and the pure DSC state where a resonance peak 
emerges for an infinitesimally small interaction strength. 
We note, however, that for the strength of the fermionic 
interaction usually taken to describe the resonance peak 
in the DSC state, a resonance peak also emerges in 
the DDW state. Moreover, we find that the resonance 
peak in the DDW state is basically dispersionless and 
confined to the vicinity of Q = (tt, tt) due to the form 
of the particle-hole continuum in the DDW state. In 
contrast, the dispersion of the resonance peak in the 
coexisting DDW and DSC state is similar to that in the 
pure DSC state, with the exception that in the vicinity 
of Q, the former exhibits a cusp. These results show 
that the detailed momentum and frequency dependence 
of the resonance peak is different in all three phases, 
the pure DDW, pure DSC and coexisting DDW and 
DSC phases. Thus, a detailed experimental study of the 
resonance peak in the underdoped cuprates permits one 
to identify the nature of the underlying phase in which 
the resonance peaks emerges. We believe, however, that 
the currently available experimental INS data do not yet 
allow an unambiguous conclusion with regards to the 
nature of the phases present in the underdoped cuprate 
superconductors. 



the particle-hole continuum in the vicinity of Q in the 
coexistence phase, as is evident from Fig. Efa). Thus 
the dispersion of the resonance peak directly reflects the 
different momentum dependence of the particle-hole con- 
tinuum in the vicinity of Q in the coexisting DDW and 
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APPENDIX A: DYNAMICAL SPIN SUSCEPTIBILITY IN THE REGIME OF COEXISTING DDW+DSC 

PHASES 



In the coexisting DDW and (i-wave superconducting phase, the susceptibility is given by 

X(q,w) = ^x«(q,w) 



(Al) 



where 



x (1) (q^) 



A k A 



k^k+q 



"k"k+q 



- r k £ k+q + WkWk+q 



/(»k +q ) - mi) 

(A2) 



x (2) (q,^) 



32 ^ ( 1 



^k^k +q + AkA k+q N 



"k"k+q 



£ k £ k+q + W k W k+tl 



l-/("k+q)-/(»k) | /("k+q) + /("k)-l ' 

'~ to-Q+ +(l -Q+ + iS 



^ + fi k+q + Q k + iS 



(A3) 



x (3) (q,^) 



^k g k+q 



AkA 



k^Ak+q 



16 



"k"k+q 



£ k £ k+q + ^k+q 



/("k + q)-/(^) 



/(»+)- /(» k+q ) \ 



W - ^ k +q + ^k + ^ W - fij + ^k+q + ^ 



(A4) 



^k ^k+q + AkAk+q 



"k"k+q 



l-/(^k +q )-/W , /M+Z^k+q)- 1 



£ k £ k+q+W k+q 



V 



V / ( £ k) 2 + ^kV( £k+q ) 2 + ^k 2 +q y' 



w + ft+ + ft k+q + iS 



" k+q - ^k + M 



(A5) 



x (5) (q,^) 



32 ? ( 1 



^k ^k+q + AkAk +q s 



^k^k+q 



£ k £ k+q + WkWk+q 



l-/OW-/(°k) , /("k)+/("k + q)-l N 



r k+q) 



u> + n k + o k+q 



W - "k+q - ft k 



•t<S 



(A6) 



x (6) (q,^) 



K ^ k +q + AkAk+q ^ 



16 V 

k 



^k ^k+q 



£ k £ k+q 



k+q 



fe) 2 +^k 2 v( £k +q) +^+ q ; 



Z("k +q )-/("k) 
w - ^k+ q + ^k + m 



(A7) 



T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999); 
M. R. Norman, D. Pines, and C. Kallin, cond-mat/0507031 
(unpublished). 

P.W. Anderson et al, J. Phys. Condens. Mat. 16, R755 
(2004); D. A. Ivanov, P. A. Lee, and X.-G. Wen, Phys. Rev. 
Lett. 84, 3958 (2000); CM. Varma, Phys. Rev. Lett 83, 
3538 (1999); V. J. Emery, S. A. Kivelson, and O. Zachar, 
Phys. Rev. B 56, 6120 (1997); L. Benfatto, S. Caprara, 
and C. Di Castro; Eur. Phys. Jour. B 17. 95 (2000); J. 
Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. Lett. 
80, 3839 (1998); J.R. Engelbrecht, A. Nazarenko, M. Ran- 
deria, and E. Dagotto, Phys. Rev. B 57, 13406 (1998); Q. 
Chen, I. Kosztin, B. Janko, and K. Levin, Phys. Rev. B 
59, 7083 (1999); S.C. Zhang, Science 275, 1089 (1997). 
S. Chakravarty, R.B. Laughlin, D.K. Morr, and C. Nayak, 
Phys. Rev. B63, 094503 (2001). 

S. Tewari, H.-Y. Kee, C. Nayak, and S. Chakravarty, Phys. 
Rev. B 64, 224516 (2001). 

S. Chakravarty, C. Nayak, S. Tewari, and X. Yang, Phys. 
Rev. Lett. 89, 277003 "(2002). 

M.V. Eremin, I. Eremin, and A. Terzi, Phys. Rev. B 66, 
104524 (2002). 

J. Rossat-Mignod, L.P. Regnault, C. Vettier, P. Bourges, 

P. Burlet, J. Bossy, Physica C 185-189, 86 (1991). 

H.A. Mook, M. Yethiraj, G. Aeppli, T.E. Mason, and T. 

Armstrong, Phys. Rev. Lett. 70, 3490 (1993). 

H. F. Fong, B. Keimer, P.W. Anderson, D. Reznik, F. Do- 

gan, and LA. Aksay, Phys. Rev. Lett. 75, 316 (1995); ibid., 

Phys. Rev. B 54, 6708 (1996); H.F. Fong et al, Nature 

(London) 398, 588 (1999); S. Pailhes, Y. Sidis, P. Bourges, 

V. Hinkov, A. Ivanov, C. Ulrich, L.P. Regnault, and B. 

Keimer, Phys. Rev. Lett. 93, 167001 (2004). 

P. Bourges, L.P. Regnault, Y. Sidis, and C. Vettier, Phys. 

Rev. B 53, 876 (1996); H. He et al, Science 295, 1045 

(2002). 

See for review P. Bourges, in "The gap Symmetry 
and Fluctuations in High Temperature Superconductors" 
edited by J. Bok, G. Deutscher, D. Pavuna and S.A. Wolf 
(Plenum Press, 1998). 

P. Dai, H.A. Mook, S.M. Hayden, G. Aeppli, T.G. Per- 
ring, R.D. Hunt, and F. Dogan, Science 284, 1344 (1999); 
P. Dai, M. Yethiraj, H.A. Mook, T.B. Lindemer, and F. 
Dogan, Phys. Rev. Lett. 77, 5425 (1996). 
H.F. Fong, B. Keimer, D.L. Milius, and LA. Aksay, Phys. 
Rev. Lett. 78, 713 (1997); H.F. Fong et al, Phys. Rev. B 
61, 14773 (2000). 

S.M. Hayden, H.A. Mook, P. Dai, T.G. Perring, and F. 
Dogan, Nature (London) 429, 531 (2004). 
C. Stock, W.J.L. Buyers, R.A. Cowley, P.S. Clegg, R. 
Coldea, CD. Frost, R. Liang, D. Peets, D. Bonn, W.N. 



Hardy, and R. J. Birgeneau Phys. Rev. B 71, 024522 
(2005); C. Stock, W.J.L. Buyers, R. Liang, D. Peets, Z. 
Tun, D. Bonn, W.N. Hardy, and R.J. Birgeneau, Phys. 
Rev. B 69, 014502 (2004). 

16 H.F. Fong et al. Phys. Rev. Lett. 75, 316 (1995); D. Z. Liu, 
Y. Zha, and K. Levin Phys. Rev. Lett. 75, 4130 (1995); A. 
J. Millis and H. Monien Phys. Rev. B 54, 16172 (1996); 
A. Abanov, and A.V. Chubukov, Phys. Rev. Lett. 83, 
1652 (1999); T. Dahm, D. Manske, and L. Tewordt, Phys. 
Rev. B 58, 12454 (1998); J. Brinckmann and P. A. Lee, 
Phys. Rev. Lett. 82, 2915 (1999); Y.-J. Kao, Q. Si, and 
K. Levin, Phys. Rev. B 61, R11898 (2000); F. Onufrieva 
and P. Pfeuty, Phys. Rev. B 65, 054515 (2002); M. Eschrig 
and M.R. Norman, Phys. Rev. Lett. 89, 277005 (2002); D. 
Manske, I. Eremin, and K. H. Bennemann, Phys. Rev. B 
63, 054517 (2001); M.R. Norman, Phys. Rev. B 61, 14751 
(2000); ibid 63, 092509 (2001); A.V. Chubukov, B. Janko 
and O. Tchernyshov, Phys. Rev. B 63, 180507(R) (2001); I. 
Sega, P. Prelovsek, and J. Bonca, Phys. Rev. B 68, 054524 
(2003). 

17 I. Eremin, D.K. Morr, A.V. Chubukov, K. Bennemann, and 
M.R. Norman, Phys. Rev. Lett. 94, 147001 (2005). 

18 L. Yin, S. Chakravarty, and P.W. Anderson, Phys. Rev. 
Lett. 78, 3559 (1997); E. Dernier and S.C. Zhang, Phys. 
Rev. Lett. 75, 4126 (1995); D.K. Morr and D. Pines, Phys. 
Rev. Lett. 81, 1086 (1998); M. Vojta and T. Ulbricht Phys. 
Rev. Lett. 93, 127002 (2004); G. S. Uhrig, K. P. Schmidt, 
and M. Griininger Phys. Rev. Lett. 93, 267003 (2004); 
G. Seibold and J. Lorenzana Phys. Rev. Lett. 94, 107006 
(2005). 

19 A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

20 E. Cappelluti and R. Zeyher, Phys. Rev. B 59, 6475 (1999); 
I. Eremin and M. Eremin, J. Supercond. 10, 459 (1997). 

21 B. Dora, K. Maki, and A. Virosztek, Mod. Phys. Lett. 18, 
327 (2004). 

22 S. Chakravarty, Ch. Nayak, and S. Tewari, Phys. Rev. B 
68, 100504(R) (2003). 

23 B. Valenzuela, E.J. Nicol, and J. P. Carbotte, Phys. Rev. 
B 71, 134503 (2005). 

24 See, for example, G.D. Mahan, Many- Particle Physics, 
(Plenum Press, New York and London, 1990). 

25 The results for Im^o were obtained by using the tetrahe- 
dron method— to evaluate Eq.@. 

26 G. Lehmann, and M. Taut, Phys. Stat. Solidi B 54, 469 
(1972). 

27 Despite the doubling of the unit cell one finds for the umk- 
lapp susceptibility x(cu Q> w ) — in the DDW state. 



